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The incoherent approximation for the determination of the vibrational density of states of glasses 
from inelastic neutron or x-ray scattering data is extended to treat the coherent scattering. The 
method is applied to new room temperature measurements of vitreous silica and germania on the 
thermal time-of-flight spectrometer IN4 at the High Flux Reactor in Grenoble. The inelastic dynamic 
structure factor at the boson peak turns out to be reasonably well described in terms of a mixture 
of rotation and translation of practically undistorted SiC>4 or GeC>4 tetrahedra. The translational 
component exceeds the expectation of the Debye model by a factor of two. A possible relation of 
this excess to the phonon shift and broadening observed in x-ray Brillouin scattering experiments 
is discussed. 

PACS numbers: 63.50.+X, 64.70.Pf 
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I. INTRODUCTION 

There is an extensive quantum mechanical treatment 
of the scattering from atoms moving in a crystal 1 . The 
regular atomic arrangement allows to solve this problem 
with an accuracy and a theoretical depth which one can- 
not hope for in the much more complex case of a disor- 
dered solid. In a glass, the translation symmetry is lost. 
In addition, the sample is in an energy landscape with 
many minima. Fortunately, the relaxational jumps be- 
tween different energy minima of glasses are only seen as 
quasielastic scattering below an energy transfer of about 
1 meV«2i£ (a frequency of 250 GHz); above that frequency, 
one can reckon with a more or less harmonic vibrational 
density of states. 

But the problem what these vibrations are is by no 
means solved. Below 1 meV, reasonably well-defined 
long- wavelength sound waves can be shown to exist (co- 
existing with relaxational or tunneling motion) , but they 
become rapidly overdamped above that frequency. This 
was first deduced* from the plateau in the thermal con- 
ductivity at low temperatures. Within the last decade, 
it has been directly observed for the longitudinal sound 
waves by x-ray Brillouin scattering 5 . Also, the density of 
states exceeds the Debye expectation markedly at the so- 
called boson peak, at a frequency where the correspond- 
ing crystals still have only well-defined long-wavelength 
sound waves. 

In this paper, we report thermal neutron time-of-flight 
measurements at room temperature on silica and ger- 
mania, data taken over a very large momentum transfer 
range with the new spectrometer IN4 at the High Flux 
Reactor at the Institut Laue-Langevin at Grenoble. The 



high quality of the data allows for an evaluation which 
goes beyond the usual incoherent approximation, extend- 
ing and improving earlier work&L&SiiS. One not only gets 
a vibrational density of states which compares favorably 
with heat capacity data, but one learns new facts on the 
details of the atomic motion, in particular in the low- 
frequency range at the boson peak. 

After this introduction, the paper describes an exten- 
sion of the incoherent approximation for the evaluation of 
coherent inelastic neutron or x-ray data in Section II. Sec- 
tion III presents the time-of-flight experiments in silica 
and germania, their multiple-scattering correction and 
their normalization to diffraction data from the litera- 
ture. Section IV applies the new method to the data. 
The discussion of the results and a short summary is 
given in Sections V and VI. 



II. EXTENDING THE INCOHERENT 
APPROXIMATION 

A. Definitions 

The following derivation is formulated in terms of the 
classical scattering law S(Q, to), where the frequency to is 
related to the energy transfer E of the scattering process 
by E = hco and Q is the momentum transfer. This classi- 
cal S(Q,u>) can be calculated from the measured double 
differential cross section 
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taking ui to be positive in energy gain of the scattered 
particle, neutron or x-ray photon. Here T is the temper- 
ature, fcj and kf are the wavevector values of incoming 
and scattered waves, respectively, N is the number of 
atoms in the beam, a is the average scattering cross sec- 
tion of the atoms and is the solid angle. Note that 
the average scattering cross section a is Q-dependent for 
x-rays; this Q-dependence is determined by the atomic 
form factors. The definition of S(Q,lo) requires a com- 
pletely isotropic glass, the usual case. It is also valid for 
glasses with more than one kind of atom, like silica and 
germania. 

The above definition makes no distinction between co- 
herent scattering (the case where the waves scattered 
from different atoms interfere) and incoherent scattering 
(the case where there is no interference). X-ray scat- 
tering is purely coherent; for neutrons, it depends on the 
nuclei of the scattering atoms. In silicon, germanium and 
oxygen atoms, the coherent scattering dominates. Inte- 
grating over all frequencies one obtains S(Q) 

/oo 
dwS(Q,uj). (2) 
-oo 

If all atoms of the sample scatter only incoherently, 
S(Q) = 1. If one has only coherent scattering, S(Q) 
shows a first sharp diffraction peak at about 1.5 A, fol- 
lowed by oscillations around 1 which die out at high Q. 
These peaks reflect the short range order of the glass on 
the atomic scale. Below the first sharp diffraction peak, 
S(Q) is due to long-range density and concentration fluc- 
tuations. Silica and germania show a pronounced first 
sharp diffraction peak. 

For the evaluation of coherent scatterers, we will need 
the definition of a hypothetical incoherent scattering 
function Si nc {Q,u). This is defined as the scattering 
function which one would have without interference be- 
tween different atoms, keeping their total cross sections. 

B. The incoherent approximation 

To determine a vibrational density of states from scat- 
tering data, one needs to take the Debye- Waller factor 
and the multiphonon scattering into account. An ele- 
gant way to do this is to use the intermediate scattering 
function 

/oo 
cos ujt S{Q,iv)duj. (3) 
-oo 

with the back transform 

i r°° 

S{Q,lu) = - cosojt S{Q,t)dt, (4) 
n Jo 

The incoherent approximation assumes that one can 
describe the scattering function, eq. Q, in terms of 
an average atom which scatters only incoherently. The 
time-dependent displacement of this average atom from 



its equilibrium position is assumed to have a Gaussian 
distribution. From the Bloch identity, one obtains the 
intermediate scattering function 

Sinc(Q,t)=e- Q2 ^\ (5) 

where 7(f) is the time-dependent mean square displace- 
ment of the average atom. 

If there are only vibrations, j(t) is determined by the 
vibrational density of states g{ui). Their relation can be 
derived from the one-phonon approximation for the in- 
elastic scattering from our classical isotropic incoherent 
scattered 

S mc (Q,cHQ 2 e- 2 ^^M ( 6) 

where e~ 2W is the Debye- Waller factor and M is the 
average atomic mass. Comparing the initial Q 2 rise and 
using the Fourier transformation, eq. one finds 

7(*) = K p^^l-cos^). (7) 

M Jo w 

Here Lu max is the upper frequency boundary of the vibra- 
tional density of states. 

The incoherent scattering is obtained from the Fourier 
transform of the intermediate scattering function in time 

Sinc{Q,u) = -\ dt cos tote-" < (t)Q2 . (8) 

7T Jo 

In this approximation, one accounts not only for the one- 
phonon scattering, but the entire inelastic scattering in- 
cluding the multiphonon terms. 

To apply the approximation to measured data, one be- 
gins by calculating a first guess to the vibrational density 
of states, assuming a Q 2 - or Q 2 e _2W '-dependence of the 
inelastic scattering. From the density of states, one gets 
7(i) from eq. J7J) and can calculate the Q-dependence of 
the inelastic scattering. With this, one can determine a 
better density of states from the data. Usually, after a 
few iterations the density of states does no longer change. 
In the case of a glass consisting of two or more elements, 
one calls this density the generalized vibrational density 
of states, to emphasize that is is not the true density of 
states, but its reflection in the scattering, weighted by 
the cross sections and masses of the atoms in the sample. 

As we will see, the incoherent approximation works 
astonishingly well even for the two coherently scattering 
glasses silica and germania. But it does not provide any 
information on the vibrational eigenvectors. One needs 
an extension of the incoherent approximation to do that. 
Such an extension will be introduced in the following sub- 
section. 

C. Oscillation function S U (Q) 

The vibrational density of states is a function of the 
frequency, not of the time. Therefore the extension of 
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the incoherent approximation must be done in the fre- 
quency domain. The vibrational eigenvectors change 
with changing frequency, so each frequency window has 
its own coherent dynamic structure factor. The interfer- 
ence between different scattering atoms does not change 
the overall scattering intensity, but leads to oscillations 
in the dependence on the momentum transfer. 

To take this into account, we define the oscillation 
function S U (Q) by 



s u (Q) = 



S(Q,u>) 
S mc (Q,u) ' 



(9) 



where Si nc (Q,uj) is understood to be the scattering func- 
tion without interference between different atoms, but 
with unchanged total cross sections, as defined in II. A. 

Like S(Q), S U) (Q) equals 1 in the incoherent case. It is 
an extension of S(Q) to describe the coherent oscillations 
of the scattering at a fixed frequency uj, both for elastic 
and inelastic scattering. Though our main interest is in 
the inelastic part, it is suitable to begin the discussion 
with the elastic part Sq(Q) for u = 0. 

As S(Q) reflects the pair correlation of the atoms in 
their instantaneous positions, the elastic oscillation func- 
tion Sq(Q) reflects the pair correlation of the atoms in 
their equilibrium positions. In a glass at low tempera- 
tures (and even at not so low temperatures), the atomic 
displacements are small compared to the interatomic dis- 
tances, so Sq(Q) « S(Q). Both functions show the first 
sharp diffraction peak at about 1.5 A followed by fur- 
ther peaks at higher momentum transfer. At still higher 
momentum transfer, the peak amplitudes diminish and 
the functions S(Q) and Sq(Q) approach the final value 
of 1. Below the first sharp diffraction peak, at small mo- 
mentum transfer, Sq(Q) mirrors the frozen density and 
concentration fluctuations of the glass, while S(Q) mir- 
rors both static and dynamic ones. 

The inelastic part of the oscillation function S UJ (Q) 
with w ^ will depend on the vibrational modes at 
the given frequency. At small momentum transfer, it 
shows the Brillouin peak, the scattering from longitudinal 
sound waves. The interference pattern at higher momen- 
tum transfer is not only due to the positional phase fac- 
tors which determine S(Q), but also to the scalar product 
of momentum transfer and atomic displacement (see eq. 
lHJl ). Thus the inelastic part of Suj(Q) contains informa- 
tion on the eigenvectors of the vibrational modes at the 
frequency ui. 

With the help of the oscillation function, the incoher- 
ent approximation, eq. (JSJ , transforms into the extended 
approximation 

S(Q,u) = S ll> {Q)- dt cos ute~^ t)Q2 . (10) 
* Jo 

The approximation allows to fit not only a density of 
states, but a frequency-dependent oscillation function as 
well. Thus, for coherent scatterers, the introduction of 
Su{Q) provides not only an enormous reduction of the 



deviation between incoherent approximation and exper- 
iment, but opens up the possibility to analyze the vi- 
brational eigenvectors^. If this analysis is successful, 
one can proceed to calculate the true vibrational density 
of states, because then one knows the weight of a given 
mode in the scattering. 

In the harmonic one-phonon approximation, one can 
express the coherent inelastic scattering 1 from a given 
normal mode in terms of its eigenvector &j , the equilib- 
rium position Yj and the scattering length bj of atom 
j, (j — 1..N) (the scattering length bj is related to the 
coherent cross section dj by Uj = 4tt6^). Within the 
one-phonon approximation, the oscillation function 
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where the angular brackets denote an average over all 
eigenmodes at the frequency u>, together with a direc- 
tional average over the momentum transfer vector Q. 
The mode normalization factor F norm is given by 
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If one is able to find the proper mode eigenvectors, one 
can go beyond the determination of a generalized vibra- 
tional density of states, because the proper normalization 
factor contains cross sections and masses and thus allows 
to determine the true vibrational density of states. This 
paper aims at such a treatment for the silica and germa- 
nia measurements described in the next section. 
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FIG. 1: IN4 spectrum (with standard corrections) from vit- 
reous silica at a scattering angle of 104 degrees, showing the 
boson peak at both sides of the elastic line. 
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III. EXPERIMENT 
A. Samples and time-of-flight experiment 

The vitreous silica sample used for the IN4 experiments 
was a commercial grade spectrosil disk with a diameter 
of 50 mm and a thickness of 4.8 mm. 

In the case of vitreous germania, appropriate amounts 
of reagent-grade GeC>2 powder (Aldrich 99.99+%) were 
melted in platinum crucibles for about 1 h at ca. 1600 C. 
The homogeneous and bubble- free melt was subsequently 
quenched in water and glassy samples of irregular shape 
were removed from the bottom of the crucible. This 
preparation technique results in completely transparent 
glasses. Several such pieces were arranged to mimic the 
disk shape of the silica sample. 

The neutron spectra were obtained on the thermal 
time-of-flight spectrometer IN4 at the High Flux Reac- 
tor of the Institute Laue-Langevin at Grenoble. The 
measurements were done in reflection geometry, with the 
sample disk plane inclined at 45 degrees to the incoming 
beam. With a wavelength of incoming neutrons of 1.53 
A, one is able to study the momentum transfer range 
from 1.5 to 7 A with an energy resolution of 1.3 meV 
FWHM. Similarly, a wavelength of incoming neutrons 
of 2.2 A allowed to study the momentum transfer range 
from 1 to 4.9 A with an energy resolution of 0.8 meV 
FWHM. 

The momentum transfer range of these measurements 
is more than a factor of two larger than the one of earlier 
investigational with cold neutrons on the time-of-flight 
spectrometer IN6. One of these earlier measurements 
(vitreous silica at 318 K) was included in the evaluation 
presented below. 

The IN4 measurements were performed at tempera- 
tures between 5 and 300 K. The following evaluation, 
however, is restricted to the 300 K data. 



self to get the spectrum of the twice-scattering processes. 
However, to treat the elastic line correctly, one has to re- 
place the elastic line by a (5-function in one of the two 
spectra which one folds. Otherwise the procedure would 
broaden the elastic line in an unphysical way. 

The question how much of the resulting twice- 
scattering spectrum one should subtract is answered by 
looking at the momentum transfer dependence in the in- 
elastic part of the spectrum. On a normalized scale, one 
subtracts a fraction s mu , chosen in such a way that the 
inelastic intensity extrapolates to zero towards zero mo- 
mentum transfer. As an example, Fig. 2 shows corrected 
(with s mu = 0.1) and uncorrected data for germania at 
300 K between 1 and 5 meV, measured with incoming 
neutrons of a wavelength of 2.2 A. 

The procedure is not exact, because the coherent in- 
elastic scattering at small momentum transfer is not ex- 
actly zero. If the velocity of the incoming neutrons ex- 
ceeds the longitudinal sound velocity of the glass (the 
kinematic condition^ for the visibility of the Brillouin 
line), one sees the Brillouin scattering. But even if it 
does not (as in the measurements reported here), there is 
still a small nonzero coherent inelastic scattering contri- 
bution. This, however, is small compared to the multiple 
scatterinssiiii^. 



C. Normalization to S(Q) 

In principle, if one does all corrections properly, one 
should be able to normalize the measurement to the vana- 
dium signal. In practice, it is easier and more accurate to 
normalize the time-of-flight data of a glass to the S(Q) of 
a diffraction measurement. For vitreous silica, there are 
two such diffraction measurements in the literatur e) 13 ! 14 . 
They agree very well with each other. For germania, 
there is only one diffraction measurement^. 

To determine S(Q) from a time-of-flight dataset, one 



B. Corrections 

The measured neutron counts were corrected for the 
empty container signal. The detector efficiency was cor- 
rected for by a measurement of a vanadium sample (an 
incoherent scatterer) in the same container. The signal 
was multiplied by ki/kf and an angle-dependent absorp- 
tion correction was calculated. These four standard cor- 
rections were done with the program INX of the Institut 
Laue-Langevin at Grenoble. Fig. 1 shows such a cor- 
rected spectrum for silica at the highest scattering angle. 

The multiple scattering correction was done assuming 
an angle-independent multiple scattering contribution. 
For silica and germania, where wide-angle scattering pre- 
dominates, this assumption is expected to hold. 

To do the correction, one first calculates an average 
spectrum, weighting each detector with the sine of its 
scattering angle. This spectrum should be folded with it- 
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FIG. 2: Inelastic signal from 1 to 5 meV for germania at 300 
K, incoming wavelength 2.2 A, with and without multiple 
scattering correction. 
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integrates the scattering over the time channels of a single 
detector or a detector group, correcting for the change 
in momentum transfer as the energy transfer changes. 
Here, we corrected with the factor Q 2 el /Q 2 , where Q e i is 
the momentum transfer at the clastic line and Q is the 
one at finite energy transfer. The procedure is not exact, 
but it has the advantage that a bad detector does not 
corrupt its neighbors. 

Fig. 3 (a) shows the result of this normalization to 
S(Q) for the two IN4 measurements on vitreous silica 
at room temperature. Also included is an earlier IN6 
measurement^, done with incoming neutrons of a wave- 
length 4.1 A at 318 K. 

The same normalization of the two IN4 room temper- 
ature measurements of vitreous germania to diffraction 
data from this substance^ is shown in Fig. 3 (b). The 
1.5 A measurement suffers from problems with the sub- 
traction of the empty container; it consisted of pure alu- 
minum with large crystalline grains. 

The comparison to Fig. 3 (a) shows that the second 
and the third peak in S(Q) shift to lower momentum 
transfer in germania. This is caused by the increased Ge- 
O distance (1.73 A) in germania^ as compared to the 
Si-0 distance of 1.6 A in silica^. The corner-connected 
Ge04-tetrahedra in germania are larger than the SiC>4 



ones in silica. Nevertheless, the position of the first sharp 
diffraction peak is more or less the same in both sub- 
stances, which shows that the packing of the tetrahedra 
must be different. 

The normalization allows to compare different mea- 
surements in the same frequency window. Fig. 4 (a) 
shows the boson peak frequency window from 2 to 6 meV 
in vitreous silica for all three sets of data. One sees pro- 
nounced oscillations around the dashed line (calculated 
from the incoherent approximation as explained in Sec- 
tion III. A). Germania shows rather similar oscillations in 
the same window in Fig. 4 (b) (its boson peak is also at 
« 4 meV) . The task of the next section is to extract the 
information content in these oscillations. 

As one goes up in energy transfer, the amplitude of 
the oscillations decreases gradually. This development 
is shown for vitreous silica in Fig. 5 (a-c). Above 40 
meV, one finds the incoherent approximation to be valid 
within the experimental error. In vitreous germania, the 
incoherent approximation is reached even faster, around 
30 meV. In both cases, there is no discernible peak shift; 
within experimental error, the peaks merely fade away 
with increasing energy transfer. 
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FIG. 3: (a) Normalization of IN4 data at 2.2 A and 1.5 A to 
the S(Q) of vitreous silicai-*— (continuous line). An earlier 
IN6 measurement- at 4.1 A is also included (b) Normaliza- 
tion of IN4 data at 2.2 A and 1.5 A to the S(Q) of vitreous 
germaniai* (continuous line). 
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FIG. 4: The inelastic dynamic structure factor in the bo- 
son peak region from (a) the three measurements of vitreous 
silica (b) the two measurements of vitreous germania. The 
dashed line is the incoherent approximation discussed in Sec- 
tion III. A, the continuous line the extension to coherent scat- 
tering in terms of the empirical five-lorentzian fit explained 
in Section III.B. 
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IV. EVALUATION 

A. Incoherent approximation 

Having data corrected for multiple scattering and nor- 
malized to S(Q), one can proceed to apply the incoher- 
ent approximation described in Section II. B. In this way, 
one determines a generalized vibrational density of states 
g{E) (we replace the frequency lo by the energy transfer 
E = Tuij in this section) without any adaptable parame- 
ter. 

In both cases, silica and germania, the procedure de- 
scribed in section II. B converged to a final generalized 
vibrational density of states after three iteration steps. 

As it turns out, the three results for g(E) of vitreous 
silica agree fairly well with each other. The boson peak in 
g(E)/E 2 is slightly lower in the 2.2 A measurement, but 
this is not due to the oscillations of the inelastic dynamic 
structure factor. Looking at Fig. 4 (a), one sees that the 
normalized intensity itself is slightly lower than in the 
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FIG. 5: The development of the inelastic dynamic structure 
factor of vitreous silica above the boson peak. The dashed line 
is the incoherent approximation discussed in Section III. A, the 
continuous line the extension to coherent scattering in terms 
of eq. tHj . (a) Energy transfer from 10 to 15 meV (b) Energy 
transfer from 20 to 30 meV (c) Energy transfer from 40 to 60 
meV. 



two other measurements. 

The same good agreement was found for the vibra- 
tional densities of states of germania calculated for the 
IN4 measurements at 1.5 and 2.2 A. 



Analysis of the inelastic dynamic structure 
factor 



As the following analysis shows, we have at present 
no perfect eigenmode fit of the dynamic structure fac- 
tor at the boson peak. Therefore we do the analysis in 
two steps. We first fit the oscillation function Sb oson (Q) 
at the boson peak in a purely empirical way with five 
lorentzians, independent of any motional model, but pro- 
viding an excellent fit. This fit form is used to investigate 
the frequency dependence of the peaks in S UJ (Q) to higher 
frequencies. The second step models the atomic motion 
at the boson peak in terms of a sum of translation and ro- 
tation of undistorted tetrahedra. In this procedure, one 
has only one fit parameter (the translational fraction of 
the motion). One gets a reasonable rather than a per- 
fect fit, but one has a motional model which allows to 
calculate the true rather than the generalized vibrational 
density of states. 

To go beyond the incoherent approximation, we start 
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FIG. 6: Inelastic dynamic structure factor Sboson(Q) in the 
boson peak region for (a) silica (continuous line sum of 0.45 
S(Q) and 0.55 tetrahedra libration, see text) (b) germania 
(continuous line sum of 0.75 S(Q) and 0.25 tetrahedra libra- 
tion) . 
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by calculating the oscillation function S LJ (Q) at the boson 
peak via 



bosoji 



(Q) 



oson 

(Q,w) 



l / tt dt cos ujte-fWQ 2 



(13) 



taking j(t) from the evaluation in terms of the incoher- 
ent approximation described in the previous subsection. 
Sboson(Qi is determined from the measured S(Q,u>) 
via 



os on, 



(14) 



where huj m i n =2 meV was chosen to exclude any contri- 
bution of the elastic line even for the 1.53 A measurement 
and tko max =6 meV is a frequency well above the two bo- 
son peak frequencies of 4 meV and 3.6 meV for silica and 
germania, respectively. Fig. 6 (a) shows the result for 
the three sets of data of silica, Fig. 6 (b) for the two IN4 
measurements of germania. 

The next step is to quantify the fading-away of the os- 
cillations with increasing frequency. For this, we need 
a functional expression for Sboson(Q)- To get it, we 
fit the data points of Fig. 6 (a) in terms of a sum 
of five lorentzians (the continuous line in Fig. 6 (a)). 
This purely empirical function oscillates around 1 in the 
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momentum transfer region of the three measurements. 
One can use it to replace the incoherent approximation 
Sboson(Q) = 1- Fig. 4 (a) shows that it gives a much 
better agreement with experiment than the incoherent 
approximation, at least in the frequency region of the 
boson peak. 

As one goes up in frequency, the oscillations of the 
measured dynamic structure factor begin to get weaker, 
until one reaches again the incoherent approximation at 
about 40 meV (see Figs. 5(a) to (c)). One can follow 
this behaviour quantitatively by fitting the inelastic in- 
tensities in subsequent frequency windows in terms of the 
oscillation function 

S U (Q) = r mcoh (E) + [1 

'incoh (E)]Sboson(Q) (15) 

where S bo son(Q) is the five-lorentzian fit of Fig. 4 and 
fincoh(E) is the energy-dependent fraction of incoherent 
scattering, going from to 1 as one goes from the boson 
peak up to higher frequencies (Fig. 7 (a)). Fig. 5 shows 
that one obtains a good fit of the measurements in this 
way. So within experimental accuracy, the peaks in the 
coherent inelastic dynamic structure factor of silica do 
not shift; they merely fade away to make room for a full 
validity of the incoherent approximation above 40 meV. 

It is interesting to note that this crossover into a valid- 
ity of the incoherent approximation occurs at the Debye 
frequency of vitreous silica. Debye's strongly over- 
simplified picture describes the whole vibrational density 
of states in terms of sound waves. The total density of 
sound waves (normalized to 1) is 



3 W 2 



(16) 
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energy transfer E (meV) 
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FIG. 7: The change of the incoherent fraction ri nco h of the 
oscillation function with increasing frequency for (a) vitreous 
silica (b) vitreous germania. 



FIG. 8: Coupled libration of five corner-connected Si04- or 
Ge04-tetrahedra®. 
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TABLE I: Debye frequencies of vitreous silica and germania 4 . 



substance 


M 


P 


Vl 


Vt 


Hu>D 




a.u. 


kg/m 3 


m/s 


m/s 


meV 


a — Si02 


20 


2200 


5800 


3800 


42 


a - GeQ 2 


34.86 


3600 


3680 


2410 


26.8 



with the Debye frequency ujd given by 



18tt 2 



J D 



M(l/vf 



(17) 



Here p is the density, vi is the longitudinal sound velocity 
and v t is the transverse sound velocity. Table I gives the 
values for vitreous silica and germania. 

In germania, the oscillations of the inelastic dynamic 
structure factor fade even more quickly with increasing 
frequency than in silica (Fig. 7 (b)). Again, the validity 
of the incoherent approximation is reached at about the 
Debye frequency. 

Having established the frequency dependence of the 
dynamic structure factor in both glasses, we proceed to 
the second step, the understanding of the oscillation func- 
tion SbosoniQ) at the boson peak in terms of a detailed 
picture of the atomic motion. 

To find such an understanding, let us first recall what 
one knows about this frequency region, both from the 
quartz crystaliS*ii and from the silica glass neutron^ 
and Hyperramani^ studies. One expects a mixture of 
long-wavelength sound waves and Si04-tetrahedra libra- 
tion. The oscillation function of the translational motion 
of long- wavelength sound wavesi^ is S(Q). The oscilla- 
tion function S ro t{Q) of the coupled tetrahedra libration 
was calculated from the motional model of Fig. 8, using 
eq. (jl 1|) for S U (Q) in the one-phonon approximation. 
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FIG. 9: Comparison of the translation oscillation function 
S(Q) (continuous line) with the oscillation function S ro t(Q) 
(dashed line) of the coupled libration of five corner-connected 
Si0 4 -tetrahedra& of Fig. 8. 



Fig. 9 compares the two oscillation functions S(Q) 
and S rot (Q) for silica. They differ mainly at the first and 
second peak. Thus a strong first peak in Sboson(Q) means 
a large sound wave fraction in the boson peak modes, a 
strong second peak a large tetrahedra libration fraction. 

Fig. 6 (a) shows the best fit (continuous line) of the 
observed Sb OSO n(Q) of silica in terms of a sum of the two 
oscillation functions. The fit is by no means perfect, but 
supplies a fraction of 0.55 of tetrahedra libration signal 
and 0.45 of translation. 

This is a surprising result. If one estimates the strength 
of the translational signal on the basis of the Debye 
sound-wave treatment, one would expect no more than a 
translational fraction of 0.2 at the boson peak^S. Simi- 
larly, the Hyperraman data— require a dominating role 
of the tetrahedra libration for their understanding. But 
obviously, the librational motion of the corner-connected 
tetrahedra is accompanied by strong translational shifts. 

A similar effect appears in germania. The fit in terms 
of a sum of the S(Q) of germania and the tetrahedra 
libration (the line in Fig. 6 (b)) gives an even larger 
fraction of 0.75 for S(Q). As we will see, the Debye 
expectation in germania is again a factor of two smaller. 
We will come back to this point in the discussion. 
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FIG. 10: Vibrational density of states (plotted as g(E)/E 2 to 
emphasize the boson peak at 4 meV), obtained using the ex- 
tended approximation for (a) the three time-of-fiight datasets 
of vitreous silica (dashed line: Debye expectation) (b) the two 
IN4 measurements of vitreous germania (dashed line: Debye 
expectation). 
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C. Vibrational density of states in the extended 
approximation 

If one knows what kind of vibrational modes one deals 
with, one can calculate the vibrational density of states 
in the extended approximation, accounting for the cross 
sections and amplitudes of the atoms participating in the 
modes. 

For vitreous silica, we take the signal at the boson peak 
to consist of 55 % of tetrahedra libration signal and 45 
% of S(Q), the result of the best fit of the data in Fig. 6 
(a). Since the tetrahedra libration is essentially oxygen 
motion, we have to correct the number of resulting modes 
dividing by the enhancement factor for a pure oxygen 
motion 



, <J M 

Jo = —irr 
a Ado 



1.49, 



(18) 



where oo is the oxygen cross section and Mo is the oxy- 
gen mass. ~5 and M are the average values of cross section 
and mass as defined in Section II. 

The increasing incoherent fraction ri nco h of Sd yn {Q) in 
eq. (I15JI is taken to be an average motion of all atoms, 
which requires no enhancement factor. 

With these assumptions, one can determine a vibra- 
tional density of states from the data using the extended 
approximation, eq. Sboson(Q) is calculated from 

the sum of 45 % of S(Q) and 55 % of the tetrahedra li- 
bration oscillation function S ro t(Q) of the model in Fig. 
8. S U (Q) is evaluated from eq. JTSJ, taking r inco h to 
follow the line in Fig. 7. 

Fig. 10 (a) shows the vibrational density of states of 
vitreous silica obtained in this way for the three sets of 
data, plotted as g(E)/E 2 to emphasize the boson peak 
at 4 meV. 

For vitreous germania, we take the signal at the boson 
peak to consist of 25 % of tetrahedra libration signal and 
75 % of S(Q). Again, we have to correct the number of 
resulting modes dividing by the enhancement factor for 
a pure oxygen motion 



f a ° M i 
Jo = — ~r-r- = 1-59. 



a Mi 



o 



(19) 



The resulting vibrational density of states g(E)/E 2 is 
shown in Fig. 10 (b). 



V. DISCUSSION 

A. Potential and limitations of the method 

Let us begin the general discussion with a disclaimer: 
The evaluation method for coherent inelastic neutron or 
x-ray scattering from glasses proposed in the present pa- 
per, which we denote by "extended approximation", is 
by no means new. Its beginnings can be traced back 
to the classical paper of Carpenter and Pelizzarrii^ (and 



even beyond that). In a less formal way, it has been 
applied earlier, not only to vitreous silica-*^, but also 
to amorphous germanium^, deuterated polybutadieneSl 
and boron trioxyde^S. Our present work merely formal- 
izes this extended approximation, introducing the con- 
cept of the oscillation function Su>(Q) and giving a recipe 
for its experimental determination. 

The oscillation function S^IQ) contains information on 
the eigenvector of the vibrational or relaxational modes 
seen at the frequency lo. The information is limited: Even 
if one knows S U (Q) with high accuracy over a large Q- 
range, one cannot determine the eigenvectors at that fre- 
quency. One can only check models of the motion against 
the measured S^iQ). In a practical sense, even such a 
check is restricted to the use of the one-phonon approx- 
imation of eq. Hll|) . because it is difficult to calculate 
the interference oscillations of the coherent multiphonon 
scattering. Fortunately, the one-phonon approximation 
holds to rather large momentum transfer at the boson 
peak, the most interesting target for these studies. 

At the boson peak, the analysis is simplified by the sci- 
entific question: one wants to know (i) whether the long- 
wavelength sound waves at the boson peak frequency still 
follow the Debye expectation (ii) whether there are ad- 
ditional modes, and if there are, what their eigenvector 
is. The first of these questions can be attacked, using the 
fac1»i£ that the oscillation function S LJ (Q) =4> S(Q) for 
long-wavelength sound waves, no matter whether they 
are transverse or longitudinal. The second requires a 
model calculation for whatever mode is expected to be 
soft in the given system. In our cases silica and germa- 
nia, one expects the coupled libration of corner-connected 
tetrahedra to be soft, because this is the soft mode of the 
phase transformation from a to j3 in crystalline quartai&. 

Note that this rather streamlined procedure has weak 
points: (1) since the sound waves interact with the 
additional modesS^Si, one needs the S U {Q) of the re- 
sulting true eigenmodes, in which translation and soft 
mode eigenvector parts have a nonzero interference term. 
Therefore their S U (Q) might be different from a simple 
sum of the two oscillation functions. (2) At the boson 
peak, measurements of incoherent scatterers2^2& reveal 
sizeable nongaussianity effects. These will tend to dis- 
tort the experimental S^iQ), determined on the basis of 
the Gaussian approximation. 

In view of these points, it is not surprising that the 
agreement between the model calculation and the data in 
Fig. 6 (a) and (b) is not perfect. Nevertheless, the com- 
parison with heat capacity data in the next subsection 
shows the reliability of the resulting vibrational density 
of states. 

The extended approximation should be particularly 
useful for the evaluation of wide-angle inelastic x- 
ray scattering data. Such measurements have been 



reportec 



but have not been evaluated on a quan- 



titative level. With the recipe given here, one could do 
a quantitative evaluation, provided one has good x-ray 
diffraction measurements. In that case, one even could 
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do only constant-Q scans (simpler to measure than the 
constant-E scans of the two referencesSL 2 ^) , relating their 
intensities by the diffraction measurement. 



B. Comparison to heat capacity 

We want to check the vibrational density of states 
obtained in the previous section from the extended ap- 
proximation against earlier results in the literature. The 
first check is against heat capacity data from silica^ and 
germaniaSi between 1 and 20 K. If one plots the heat 
capacity C p as C p /T 3 , one gets a close correspondence 
to the plot of g(E)/E 2 , showing the boson peak of silica 
and germania at about 10 K. 

In this comparison, the mode eigenvector assignment 
at the boson peak is checked, because the assumed frac- 
tion of tetrahedra libration provides the correction to the 
number of vibrational modes. 

Fig. 11 (a) compares measured heat capacity data 6 
of vitreous silica to the result of the evaluation of the 
neutron data in the extended approximation, described 
in the preceding section. We chose the results from 
the measurement at 4.1 A, because the heat capacity 
measurements^ stem from the same sample. As it turns 
out, the correction of the enhancement factor of eq. I|18|) 
is essential to get good agreement. So the comparison 
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FIG. 11: Comparison of the heat capacity calculated from 
the neutron data with and without correction for the en- 
hancement factor of the tetrahedra librational modes (the one 
without correction is close to the result from the incoherent 
approximation) to measured data (a) of vitreous silica^ (b) of 
vitreous germaniap-. 



corroborates the assignment of 55 % pure oxygen signal 
and 45 % translational signal at the boson peak in vitre- 
ous silica. 

Similarly, Fig. 11 (b) shows good agreement between 
measured 2 - 9 and calculated heat capacity data in vitreous 
germania. Here, we compare to the 2.2 A measurement, 
because it has the better resolution. Again, the correc- 
tion improves the agreement. However, in this second 
case the correction is smaller, because according to the 
fit of the measured dynamic structure factor of Fig. 6 
(b) we have only 25 % tetrahedra libration at the boson 
peak. Again, this conclusion is supported by the heat 
capacity data. 



C. Comparison to simulation 

There is a large number of molecular-dynamics sim- 
ulations of vitreous silica in the literature 3 ^* 3 ^ 3 ^, most 
of them done with effective classical force fields like the 
BKS potential 33 -. The BKS potential reproduces the 
measured S(Q), but a recent comparison to an ab ini- 
tio calculation 3 ^ shows that it fails not only to reproduce 
the vibrational density of states, but also the mode eigen- 
vectors. 

Fig. 12 compares the vibrational density of states de- 
termined from the neutron measurement at 4.1 A in the 
extended approximation with the ab initio calculation 34 
and with an earlier neutron determination with very 
short wavelength at the spallation source at Argonne 7 , 
Note that the two neutron determinations are com- 
plementary: the spallation source measurement suffers 
from an overcorrection at low frequencies, but provides 
a true picture of the high-frequency density of states. 
In contrast, the cold neutron measurement is unable to 
measure above 100 meV, but provides good results at 
low frequency, as shown by the comparison to the low- 
temperature heat capacity data in the previous subsec- 
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FIG. 12: Vibrational density of states of vitreous silica from 
an ab initio simulation 34 , from earlier neutron work 7 with 0.6 
A and from the present evaluation of the 4.1 A data. 
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tion. Between 20 and 100 meV, the two sets of data agree 
reasonably well with each other. 

Note that the ab initio simulation is unable to describe 
the boson peak region below 10 meV. Above that re- 
gion, there is impressive agreement between simulation 
and neutron experiment. But below 10 meV, the sim- 
ulated modes fail to come down to the low frequencies 
where both neutrons and heat capacity report them to 
exist. The reason might be either the small size of the 
simulation cell or its poor equilibration, natural disad- 
vantages of an ab initio simulation. 

It is interesting to compare our results to a 
simulation^* of dynamic structure factors of a different 
system, soft spheres interacting with a repulsive 1 /r 6 po- 
tential, a model appropriate for metallic glasses. In this 
case, the oscillation function at the boson peak shows 
only the peaks of S(Q). The additional modes seem to 
be a motion along a chain of nearest atomic neighbors, 
with an oscillation function which resembles S(Q). This 
is obviously different in silica and germania, where the 
resonant boson peak modes seem to have a distinct tetra- 
hedra rotation component. 



D. Sound waves at and above the boson peak 

The preceding two subsections demonstrated the abil- 
ity of the extended approximation to determine a reli- 
able vibrational density of states, in particular in the fre- 
quency region of the boson peak. This good agreement 
supports the interpretation of the boson peak modes as a 
mixture between long-wavelength sound waves and soft 
modes, in the cases of silica and germania librational 
modes of corner-connected tetrahedra. 

But the sound-wave fraction determined from the dy- 
namic structure factor is decidedly higher than the ex- 
pectation on the basis of the sound-wave Debye model. 
Looking at the boson peak region between 2 and 6 meV 
in Fig. 10 (a) and (b) and the dashed line of the Debye 
expectation, one would expect no more than 20 % sound 
waves in silica and no more than 40 % in germania. The 
fit of Sboson(Q) gives about twice as much in both sub- 
stances. This finding goes beyond earlier cold-neutron 
work<§i*, which could not quantify the sound wave frac- 
tion at the boson peak. 

The effect is apparently not limited to silica and ger- 
mania. Earlier decompositions of S(Q,uj) at the boson 
peak window into a Q 2 - and a Q 2 S(Q)-p&rt in deuterated 
polybutadiene2i and in vitreous BiO-^ also observed a 
larger Q 2 S(Q)-p&rt than expected on the basis of the 
Debye model. 

If the Q 2 S(Q) part is only due to sound waves, this im- 
plies a downward shift in frequency of the sound wave in- 
tensity above the boson peak. For the longitudinal sound 
waves, such a downward shift, together with a strong 
broadening, is in fact observed experimentally in x-ray 
Brillouin scattering^ 5 .. The broadening and the shift in- 
crease with the square of the phonon wavevector, i.e. 



with the square of the nominal frequency of the phonon. 

There is general agreement that the downward shift is 
not due to a dispersion of the sound velocity, but rather 
to sound wave scattering. One finds no evidence for any 
dispersion of either the longitudinal or the transverse 
sound velocity in tunneling diode 3 ^ or ballistic phonon- 
pulse^i experiments in vitreous silica and other glasses. 

Also, there seems to be general agreement that the 
broadening of the sound waves observed in x-ray Bril- 
louin scattering is not due to a true anharmonic damp- 
ing of the sound waves, but rather to a deviation of the 
true eigenvectors from a purely sinusoidal displacement 
in spaced. 

On the other hand, there is a hot debate on the 
proper description of the x-ray Brillouin data, whether 
one should use a damped harmonic oscillatosS (DHO) 
or whether one should take another form£& which brings 
no intensity down to the frequency zero. But both forms 
shift the intensity down to lower frequencies, consistent 
with our observation of a heightened Q 2 5(Q)-component 
at the boson peak (heightened with respect to undamped 
Debye phonons). 

The low-temperature plateau in the thermal 
conductivity* and the phonon-pulse experiments 37 
suggest that the transverse phonons have essentially 
the same fate as the one observed for the longitudinal 
phonons in x-ray Brillouin scattering 35 - 38 . 

To test these ideas, we submitted a Debye density of 
states with a Debye frequency corresponding to an energy 
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FIG. 13: Experimental density of states, plotted as g(E)/E 2 
as in Fig. 10, compared to the simple Debye expectation 
(dashed line) and to Debye sound waves with DHO damping 
as explained in the text (continuous line) for (a) vitreous silica 
(b) vitreous germania. 
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transfer of 42 meV (the value for silica) to a DHO damp- 
ing with the parameters of Benassi et a& The trans- 
verse phonons were assumed to have the same damping 
as the longitudinal ones at the same frequency, increas- 
ing proportional to the square of the phonon wavevector. 
The resulting effective Debye density of states is shown 
as the continuous line in Fig. 13 (a). To do the same 
calculation for germania, shown in Fig. 13 (b), we used 
recent neutron Brillouin measurements 39 . According to 
them, the strong damping condition T = viq (q phonon 
wavevector) is reached at the energy transfer 19.3 meV, a 
bit less than three quarters of the Debye frequency. From 
Benassi et aS^, the same condition is reached in vitreous 
silica already at 12 meV, at about only one quarter of 
the Debye frequency. 

Fig. 13 shows that the scattering of the sound waves 
is strong enough to have a non-negligible effect on their 
spectral appearance in a scattering experiment. But it 
also shows that one cannot explain the boson peak in 
terms of the sound wave scattering alone. This does not 
depend on the choice of the DHO. If one takes the al- 
ternative proposed by the Montpellier group2£, one does 
indeed get a peak, because this alternative shifts the in- 
tensity to a finite frequency, not to the frequency zero. 
But the intensity of the peak remains too small. What 
one really needs is a mechanism which raises the total 
vibrational mean-square displacement by nearly a factor 
of two over the Debye expectation, a mechanism which 
brings vibrations from higher frequencies down to the bo- 
son peak. In vitreous silica and vitreous germania, one 
has the additional advantage that one can identify these 
additional vibrations by their dynamic structure factor. 



VI. SUMMARY 

We introduce a formal treatment of the coherent inelas- 
tic neutron or x-ray wide-angle scattering from glasses, 
which takes the interference oscillations explicitly into 
account. This "extended approximation" is an extension 
of the incoherent approximation. It allows to fit a newly 
defined "oscillation function" S UJ (Q) at each frequency 
u, thus supplying information on the vibrational eigen- 
vectors. The method should be particularly useful for 
the quantitative evaluation of wide-angle inelastic x-ray 
scattering measurements. 

The application of the method to new room- 
temperature thermal neutron time-of-flight measure- 
ments of silica and germania not only provides a vibra- 
tional density of states in excellent agreement with heat 
capacity and simulation data, but also allows to quan- 
tify for the first time the ratio of tetrahedra rotation 
and tetrahedra translation at the boson peak. One finds 
about twice as much translation as in a simple Debye 
expectation. This excess is probably connected to the 
heavy damping of the sound waves at frequencies above 
the boson peak, observed by x-ray Brillouin scattering. 
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